An Interacting Galaxy System Along a Filament in a Void 
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ABSTRACT 

Cosmological voids provide a unique environment for the study of galaxy formation and 
evolution. The galaxy population in their interior have significantly different properties than 
average field galaxies. As part of our Void Galaxy Survey (VGS), we have found a system of 
three interacting galaxies (VGS_31) inside a large void. VGS_31 is a small elongated group 
whose members are embedded in a common HI envelope. The HI picture suggests a filamentary 
structure with accretion of intergalactic cold gas from the filament onto the galaxies. We present 
deep optical and narrow band H Q data, optical spectroscopy, near-UV and far-UV GALEX and 
CO(l-O) data. We find that one of the galaxies, a Markarian object, has a ring-like structure and 
a tail evident both in optical and HI. While all three galaxies form stars in their central parts, the 
tail and the ring of the Markarian object are devoid of star formation. We discuss these findings 
in terms of a gravitational interaction and ongoing growth of galaxies out of a filament. VGS_31 
is one of the first observed examples of a filamentary structure in a void. It is an important 
prototype for understanding the formation of substructure in a void. This system also shows 
that the galaxy evolution in voids can be as dynamic as in high density environments. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: kinematics and dynamics — 
galaxies: structure — large-scale structure of universe — radio lines: galaxies 



1. Introduction 

Voids are vast regions occupying most of the 
volume in the universe with sizes in the range of 20 
- 50/i _1 Mpc, usually roundish in shape and largely 
devoid of galaxies (see van de Weygaert & Platen 
2011, for a recent review). In the large scale struc- 
ture of the universe we observe today, the most 
striking features along with the voids are clusters 
and filaments. In this picture, galaxies are dis- 
tributed in a filament-dominated web-like struc- 
ture. Filaments connect clusters to each other 
and, while tenuous, act like bridges (Zcl'dovich 
1970; Shandarin & Zeldovich 1989; Bond ct al. 
1996; Colbcrg ct al. 2005a; Aragon-Calvo et al. 
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2010b). From recent redshift surveys like the 
second Center for Astrophysics Redshift Survey 
(Davis ct al. 1983), the 2dF Galaxy Redshift Sur- 
vey (Collcss ct al. 2003), the Sloan Digital Sky Sur- 
vey (SDSS) (York ct al. 2000) and 2MASS redshift 
survey (Huchra et al. 2012) we see how filaments, 
bridges and sheet-like structures form substruc- 
tures and surround the undcrdcnsc regions. 

Notwithstanding the very low density of the 
void regions, we do find a dilute population of 
galaxies in their interior. These void galaxies ap- 
pear to have significantly different properties than 
average field galaxies. Previous studies based on 
redshift surveys, have shown that the void galax- 
ies are in general small, star forming blue galax- 
ies. They have a later morphological type and 
have higher specific star formation rates than the 
galaxies in average density environment. Largely 
unaffected by the complexities and processes mod- 
ifying galaxies in high-density environments, the 
galaxies living in the isolated void regions are ex- 
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pected to have retained vital clues to their for- 
mation and evolution. It has made the study of 
the relation between void galaxies and their sur- 
roundings an important aspect of the recent in- 
terest in environmental influences on galaxy for- 
mation (Szomoru et al. 1996; Kuhn et al. 1997; 
Popescu et al. 1997; Karachentseva et al. 1999; 
Grogin & Geller 1999, 2000; Peebles 2001; Hoyle & 
Vogeley 2002; Hoyle et al. 2002; Rojas et al. 2004, 
2005; Tikhonov & Karachentsev 2006; Patiri ct al. 
2006a,b; Ceccarelli et al. 2006; Wegner & Grogin 
2008; Stanonik ct al. 2009; Krcckcl ct al. 2011, 
2012). 

Void galaxies may be the rare probes of the faint 
and tenuous substructure that hierarchical struc- 
ture formation theories predict to exist in voids 
(Dubinski et al. 1993; van de Weygaert & van 
Kampen 1993; Sahni ct al. 1994; Shcth & van de 
Weygaert 2004; Furlanetto & Piran 2006; Einasto 
et al. 2011; Aragon-Calvo & Szalay 2013). Cos- 
mological simulations show how voids are hlled by 
low-density dark matter filaments, creating a net- 
work of tenuous substructures within their interior 
(van dc Weygaert & van Kampen 1993; Gottlober 
et al. 2003; Colbcrg et al. 2005b; Springcl et al. 
2006). This may indicate that the galaxies residing 
in voids are formed along these dark matter fila- 
ments, given that the simulations reveal that dark 
matter haloes are forming along them. In fact, 
some earlier observational studies have found indi- 
cations for such filamentary substructure in voids 
(Szomoru et al. 1996; Zitrin & Brosch 2008). For 
example, the latter argue that the dwarf galax- 
ies in their local galaxy sample are located on a 
dark matter filament that itself is located in a low 
galaxy density region and is accreting intergalactic 
cold gas onto the filament. 

In this study, we present the most outstanding 
example of such a filamentary void galaxy config- 
uration, VGS_31, which was found within the con- 
text of the "Void Galaxy Survey" (VGS) (Kreckel 
ct al. 2011, 2012). 

We are conducting a multiwavclength survey of 
60 void galaxies, called "The Void Galaxy Survey" 
(VGS) (Krcckcl et al. 2011). Galaxies in the VGS 
have been selected from the Sloan Digital Sky Sur- 
vey Data Release 7 (SDSS DR7) using purely ge- 
ometric and topological techniques. The sample 
was selected on the basis of galaxy density maps 
produced by the Delaunay Tessellation Field Esti- 



mator (DTFE,Schaap & van de Weygaert (2000); 
van de Weygaert & Schaap (2009)) and the sub- 
sequent application of the Watershed Void Finder 
(WVF, Platen et al. (2007); for a more general ap- 
plication of the watershed transform to the struc- 
tural analysis of the cosmic web see Aragon-Calvo 
et al. (2010a). The combination of DTFE maps 
with WVF detected voids allow us to identify the 
void galaxies from the deepest interior regions of 
identified voids in the SDSS redshift survey. The 
goal of this survey is to study the galaxy prop- 
erties in the most underdense and most pristine 
environments in the universe where the evolution 
of galaxies is expected to progress more slowly and 
relatively undisturbed. 

Our geometrically selected sample consist of 
small galaxies, with stellar mass less than 3 x 10 10 
Mq. Most of these are small, blue star forming 
disk galaxies and many of them have companions 
and extended HI disks, which are often morpholog- 
ically and kinematically disturbed (Kreckel et al. 
2011, 2012). 

We have found a system of three linearly 
aligned galaxies, VGS_31, in a void as part of 
the VGS (Figure 1 and 2). A remarkable feature 
of VGS_31 is that the whole system is embed- 
ded in a common HI envelope (Figure 3) and the 
three galaxies are at almost the same velocity 
(Figure 4). The fact that there is a small veloc- 
ity gradient throughout whole HI envelope, from 
the far east of VGS_31b to far west to VGS_31c, 
suggests that this is a filament in which the three 
galaxies are embedded. 

The system exhibits strong signs of interactions 
(Figure 1) and star forming activity with signs 
of starbursts (Figure 5 and 6). VGS_31 consists 
of a central galaxy VGSJJla and two compan- 
ions; VGS_31b and VGSJBlc. VGS_31a is opti- 
cally slightly disturbed. It has a bar like struc- 
ture and all the star formation activity is concen- 
trated there. VGS_31b is a Markarian galaxy (Mrk 
1477). Markarian galaxies are known to have UV 
continuum excess in their spectra. They have rela- 
tively high star formation rates and many of them 
contain AGNs and starburst nuclei. They span a 
wide luminosity range between -23 and -13 mag 
and have a broad range of morphologies. One in- 
teresting statistic is that an unusual number of 
Markarian galaxies occur in tight pairs or inter- 
acting systems. In this context, VGS_31b is a typ- 
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Fig. 1. — R band negative image of VGS_31. From left to right: VGSSlb: The most remarkable member of 
the system, a Markarian galaxy, has a tail and a ring. Close up images show the inner structures such as the 
bar. VGSSla: A disk galaxy with a bar structure. VGSSlc: Smallest member of the system is optically 
undisturbed. The black bar on the top- right corner represents 10 kpc (~23"). 
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Fig. 2. — The large scale structure distribution of galaxies within 13 Mpc around the VGS_31 system (in 
black) for different projections in right ascension, declination, and redshift space. The line indicates the 
orientation of the system on the plane of the sky. Surrounding galaxies are color coded by g-r color, to be 
red if g-r> 0.6 and blue if g-r < 0.6. The symbol size indicates luminosity, with larger symbols if M r < -18 
and smaller symbols if M r >-18. 
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ical example of Markarian galaxy. VGS_31b is a 
starburst galaxy and has enhanced star formation 
at center of the disk mostly concentrated in the 
bar. It has a tail, visible both in optical and in HI 
and a ring like structure around the disk. There is 
no sign for ongoing star formation activity neither 
in the tail nor in the ring (Figure 5). VGS_31c is 
significantly smaller than the other two galaxies 
and forms stars in its central part as well. 

At first glance VGS_31 looks like a normal inter- 
acting system. However the observed properties 
described above indicate a more complex picture 
as we may be witnessing the growth of structure 
along a filament within a large cosmological void. 
In fact, in an accompanying paper by Rieder et al. 
2013 (subm.) we have explored the dynamical evo- 
lution of the growth of systems resembling VGS_31 
in voids, within the context of the ACDM cosmol- 
ogy. In this study we analysed the high-resolution 
ACDM simulation Cosmo-Grid (Portegies Zwart 
et al. 2010) to see how dark matter halo systems 
similar in mass, size and environment to VGSJ31 
came to be. We found eight systems as suitable 
candidates for harbouring a VGS_31 like system 
and then investigated their merger histories. We 
found that while VGS31-like systems have a large 
variation in formation time, the environment in 
which they are embedded evolved very similarly. 
It seems to suggest that we may be witnessing the 
assembly of a filament in a void by bringing to- 
gether several smaller filamentary structures, each 
populated by individual haloes. It conjures up the 
interesting question whether the galaxies in the 
VGS_31 configuration were formed at the same lo- 
cation or whether they each originate from a dif- 
ferent location. 

In order to study the current properties and 
assembly history of the VGSJ51 system in detail, 
we have obtained a multiwavelength data set. To 
investigate the low surface brightness features, we 
performed deep B and R band imaging, and to de- 
rive star formation properties and star formation 
history we have used deep H Q and UV imaging. 
Gas morphology, kinematics and molecular hydro- 
gen content have been investigated using 21cm HI 
and CO observations. 

This paper is organized as follows: Section 
2 describes the observations, data analysis and 
the derivation of star formation properties for 
VGS_31. Section 3 gives the main observational 



results such as star formation rate (SFR) prop- 
erties, gas content and the morphology. In Sec- 
tion 4 we discuss the results. Finally, in Section 5 
we speculate on the nature and the importance of 
VGS_31. 

2. Observations and Data Reduction 

VGS_31 has been observed in several wave- 
lengths with various telescopes between 2009 and 
2012. B & R band imaging has been gathered with 
the Isaac Newton Telescope (INT) at La Palma us- 
ing the Wide Field Camera (WFC). Long slit spec- 
tra were obtained at the William Herschcl Tele- 
scope (WHT) using the Intermediate dispersion 
Spectrograph and Imaging System (ISIS). Narrow 
band H Q imaging has been done using the Hiltncr 
Telescope at the Michigan-Dartmouth-MIT Ob- 
servatory (MDM). Near UV (NUV) and far UV 
(FUV) images have been taken from the Nearby 
Galaxy Atlas (NGA) of GALEX. Radio observa- 
tions in the 21cm HI line were performed with the 
Westerbork Synthesis Radio Telescope (WSRT) 
and CO (1-0) emission spectra have been obtained 
with IRAM. A summary of the observational stud- 
ies is given in Table 1. Parameters, quantitative 
results such as SFRs, HI and molecular hydrogen 
and stellar masses are presented in Table 2, Ta- 
ble 3, Table 4 and Table 5, respectively. 

2.1. HI imaging 

We have imaged the HI in VGS_31 as part of 
the VGS project, the details of which are described 
in Kreckel et al. (2012). Observations were done 
with the Westerbork Synthesis Radio Telescope 
(WSRT) in the maxi-short configuration provid- 
ing an angular resolution of 19" x 32". The 36' 
full width half maximum of the WSRT primary 
beam is sufficient to cover the entire VGS_31 sys- 



Table 1: Observing summary for VGSJ31 
Telescope Observation details 

INT B k R band imaging 
MDM 2.4m H a imaging 

WHT Long slit spectrum 

WSRT 21cm (HI) imaging 

IRAM CO(1-0) observations 

GALEX NUV & FUV imaging 
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R.A, (2000.0) 

Fig. 3. — The HI intensity map of VGS31. HI column density intensity contours start at 1.6 xl0 19 c77i~ 2 
and increase with 4 x 10 19 cm -2 increments. Note that members of VGS_31 are aligned along a HI filament 
and appear to be embedded in a common HI envelope. 




r.a. (2000.0) relative position E — > W (arcmin) 

Fig. 4. — Position velocity (PV) diagram of VGS_31. Left: Slice position. The black line on the total HI 
map indicates the position of the slice used to create the position velocity (PV) diagram. Black cross is 
the the zero point (RA: 13:16:10). Right: PV diagram. The PV diagram is created by taking a slice along 
the whole HI structure, from the beginning of the tail to the end of VGS_31c, connecting each points. Zero 
corresponds to the red dot in the HI image on the left. This plot shows the velocity structure of the tail and 
its connection to gas around VGS_31b. Also, it shows the gas between VGS_31b and VGS_31a and VGS_31a 
and VGS_31c, respectively. Notice the change in velocity width between VGS_31b and VGS_31a. 
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Table 2 
Parameters for VGS_31 



Name SDSS ID ra 

(J2000) 


dec z m r 
(J2000) 




g-r 


TUB 


M B 


5 


VGS_31a J131606.19+413004.2 13:16:06.19 - 
VGS_31b J131614.69+412940.0 13:16:14.69 - 
VGS_31c J131559. 18+412955.9 13:15:59.18 - 


f41:30:04.25 0.021 14.75 
f41:29:40.05 0.021 14.38 
f41:29:55.96 0.021 16.78 


-20.01 
-20.38 
-17.98 


0.32 
0.50 
0.21 


14.633 
14.632 
16.735 


-20.194 
-20.193 
-18.093 


-0.64 
-0.64 
-0.64 



Note. — Column 1: Galaxy names. Column 2: SDSS IDs of the galaxies. Column 3 & 4: Right 
ascensions and declinations. Column 5: Spectrophotometric redshifts. Column 6 & 7: Apparent and 
absolute r magnitudes drawn from the SDSS DR7 and corrected for galactic extinction. Column 8: 
g — r colors drawn from the apparent model magnitudes as measured by the SDSS DR7. Column 9 & 
10: Apparent and absolute B magnitudes derived from INT B band imaging and corrected for galactic 
extinction. Column 11: <5 gives the filtered density contrast at Ry = lh Mpc as described in Kreckel 
ct al. (2011). 



tern in a single pointing. We observed 512 chan- 
nels within a total bandwidth of 10 MHz, giving 
a Hanning smoothed velocity resolution of 8.6 km 
s . Images for this paper were made with natural 
weighting to maximize sensitivity and CLEANed 
down to 0.5 mjy beam^ 1 (~1 a), reaching column 
density sensitivities of 2 x 10 19 cm~ 2 . 

2.2. Broad band imaging and photometric 
calibration 

We used the Wide Field Camera (WFC) at the 
2.1m INT for imaging in both B and R bands 
with Harris B and R filters. Total exposure times 
were 2400 second for B and 1800 seconds for the 
R band, spread over 4 exposures for the purpose 
of dithering and facilitating cosmic ray detection. 
Standard star fields were observed each night for 
the photometric calibration. Flat field exposures 
were taken at twilight at the beginning and/or end 
of each night. The data have been reduced us- 
ing the standard IRAF 1 and Photom Data Reduc- 
tion Package (STARLINK) 2 procedures for CCD 
imaging. All the optical images were trimmed and 
overscanncd followed by bias subtraction and flat 
fielding. After that all images from each filter were 
aligned and median combined. The same proce- 
dures have been followed for the standard star ob- 
servations. 



1 http:/ /iraf.noao.edu/ 
2 http: / /star-www. rl.ac.uk/ 



2.3. H a imaging and photometric calibra- 
tion 

H Q imaging has been done with the Echelle 
CCD in direct mode at the 2.4 m Hiltner Tele- 
scope. A redshifted H a filter centered at 6693 A 
has been used. To provide a measure of the contin- 
uum, R band imaging has been performed for each 
object. The total integration times for H a and for 
the continuum have been spread over 3 exposures 
for the purpose of dithering and for facilitating 
cosmic ray detection. Spectrophotomeric calibra- 
tion stars have been chosen either from Massey 
ct al. (1988) or Okc (1990). 

After performing the standard CCD reduction 
steps described in section 2.1, each combined H Q 
image has been divided by 600 and R band con- 
tinuum image by 120 in order to normalize them 
to 1 second. The mean has been calculated for an 
empty region in each image and the ratio of these 
means has been taken as the scaling factor for scal- 
ing the continuum image before subtraction from 
the H Q image. Photometric calibration of the final 
H a images has been performed following the steps 
described in Gavazzi et al. (2006) and the refer- 
ences therein. Corrections for the atmospheric ex- 
tinction and the airmass have been performed in 
the standard way, where each spectrophotometric 
calibration star observation has been fitted using 
airmass and instrumental magnitudes to get the 
atmospheric extinction coefficient. 

Contribution of the [N II] line to the observed 
flux has been estimated using the expression from 
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Kennicutt et al. (2008) and has been subtracted 
from the H Q flux. 

Correction for the foreground extinction has 
been derived from Balmer decrements, following 
the recipes in Calzetti et al. (2000) and Dommgucz 
et al. (2012) .the H^/H^ ratio ratios have been ob- 
tained from the MPA-JHU catalog for the SDSS 
DR7 3 measured through 3" fibers. We have cal- 
culated E(B-V) from these H a /H^ ratios, using 
the reddening curve from Calzetti et al. (2000) 
to obtain the corresponding extinction. By us- 
ing Balmer decrements we correct the foreground 
extinction along the entire line of sight including 
galactic extinction. Here an important point is 
that this correction assumes that the H Q /H,3 ratio 
is the same throughout the entire H Q emission re- 
gion. We checked this assumption using our WHT 
long slit spectra which measure the H a /Hp ratios 
along slits which covers the emission along both 
the minor and the major optical axis. In Table 3, 
we give both the corrected and uncorrected H a 
fluxes. 

2.3.1. SFR fromH a imaging 

Star formation rates (SFR) from the H a emis- 
sion have been calculated following the conversion 
from Kennicutt et al. (2009) where we used the 
SFR conversion factor based on a "Kroupa" ini- 
tial mass function (Kroupa & Weidner 2003): 

SFR(H a ) [Moyr" 1 ] = 5.4 x 10~ 42 x L(H a ), 

where L{H a ) is the luminosity, calculated as 

L(H a ) [erg s" 1 ] = 4 tt L> 2 (3.086 x 10 24 ) 2 I(H a ), 

where D is the distance to the galaxy in Mpc and 
I(H a ) is foreground extinction corrected flux. 

A detailed description of the H a photometry, 
foreground extinction correction and star forma- 
tion derivation will be given in a separate paper. 

2.4. Spectroscopy 

We used the Intermediate dispersion Spectro- 
graph and Imaging System (ISIS) at the 4.2 m 
WHT at La Palma to take high resolution spec- 
tra. The R1200R grating in the red arm and the 



3 The MPA-JHU catalog is publicaly available 
and may be downloaded at http://www.mpa- 
garching.mpg.dc/SDSS/DR7/archive 



R600B grating in the blue arm have been used, 
giving resolutions of respectively 0.026 nm/pixel 
over 620-720 nm and 0.045 nm/pixel over 360-540 
ran. CuNcAr lamp exposures were taken for wave- 
length calibrations. All the reduction has been 
performed using IRAF. In addition to the basic 
steps explained above, the illumination function 
along the slit has been determined via sky and 
lamp flats. Flux calibration has been done using 
spectrophotometric standard stars. 

Figure 7 shows the H Q emission along the 
optical major and minor axis of VGS_31a and 
VGS_31b. In this figure, we compare the H Q 
emission line profiles with the H Q images of the 
corresponding regions on the galaxies. 

2.5. GALEX UV observations 

GALEX NUV and FUV data have been ob- 
tained from the Nearby Galaxy Atlas (NGA), 
taken in 2004 with an exposure time of 3754 sec- 
onds in the NUV and 3002 seconds in the FUV. 
VGS_31 falls as a background object in the corre- 
sponding science frame. The data was calibrated 
on the basis of the GALEX pipelines. 

2.5.1. SFR from NUV 

The SFR has been calculated from the GALEX 
NUV and FUV luminosities and corrected for in- 
ternal dust attenuation following the method out- 
lined in Schiminovich et al. (2010). 

SFR = L uvfuv(y° un 9)W° AAuv 
Vuv 

where Ljjy is the luminosity in erg s~ 1 Hz~ 1 , 
fuv (young) is the fraction of light that originates 
in young stellar populations, rjuy is the conversion 
factor between UV luminosity and recent-past- av- 
eraged star formation rate and A\jy is the dust 
attenuation. 

To determine the values of Apuv and Anjjv 
we assume fuv{y° un 9) — 1 an d r\uv — 10 28 ' 165 
(Schiminovich et al. 2010). 

2.6. CO (1-0) emission line observations 

CO(1-0) line observations have been carried out 
with the IRAM 30-m telescope at Pico Veleta, 
Spain in 2011. We used the Eight Mixer Receiver 
(EMIR) to observe simultaneously the CO(1-0) 



8 



B Band m 

m 9 * 

10 kpc 


Halpha 

10 kpc 


GALEX FUV 

10 kpc 


GALEX NUV 

i 

10 kpc 



Fig. 5. — VGS_31 in four wavelength regimes. From top-left to bottom-right: a) INT B band image, b) MDM H tt 
continuum-subtracted image c) GALEX Far UV image, d) GALEX Near UV image. All the images are shown at the 
same physical scale. The white contours on each image indicate the corresponding H a emission regions. Note that 
there is no Hq, emission in the ring or in the tail of VGS_31b. Also note that the H Q emission is confined to the central 
parts of the galaxies, mostly in the bar structures of VGS_31a and VGS_31b. 
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Fig. 6. — BPT diagram. The BPT diagram shows the ratio of emission line fluxes of [O IIIJ/H^ to [N 
II] /H Q . The comparison sample of emission-line galaxies has been constracted from the emission- line galaxy 
sample of Terlevich et al. (1991), galaxies defined as starbursts in SIMBAD and Markarian galaxies defined 
as starbursts in Coziol (2003). This diagram is adapted from Raimann et al. (2000). 
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(rest frequency, 115.271 GHz) and the C0(2-l) 
emission line (rest frequency, 222.8118 GHz) with 
a resolution of 5km/s. The full width half maxi- 
mum (FWHM) is, respectively, -22" and -11" at 
the two frequencies. The Wideband Line Multi- 
ple Autocorrelator (WILMA) and the Fast Fourier 
Transform Spectrometer (FTS) were used as back- 
ends. WILMA and FTS cover a channel width of 
4 GHz with 2MHz and 195 KHz resolution, re- 
spectively. The observations were carried out in 
wobbler switching mode with a frequency of 1Hz 
and a throw of 120 ". The data were reduced with 
the CLASS software. 

We have detected CO(l-O) emission from 
VGS_31b and the profile is shown in Figure 7. 
VGS_31b was not detected in CO(2-l). The de- 
tection limit for CO(2-l) is ~ 1 K km s" 1 . This is 
consistent with a normal CO(1-0)/CO(2-1) ratio, 
provided the CO(2-l) emission is distributed over 
the CO(l-O) beam area. 

The data for VGS_31a was not usable because 
of an error in the focus setting. VGS_31c has not 
been observed. 

2.6.1. CO Luminosity and Mass 

IRAM spectra are expressed in terms of an- 
tenna temperature (T^). In order to convert it to 
a flux density we adopted the expressions follow- 
ing Costagliola et al. (2011) and Saintonge et al. 
(2011): 

F 

ScoUy] = T* A [K] x 3.906 x -^l, 

' VA 

where F e ff is the forward efficiency and r\A is the 
antenna efficiency. The total flux is then derived 
by integrating the observed CO profile over veloc- 
ity. 

From the total flux we have calculated the lu- 
minosity as follows: 

L' co = 3.25 x 10 7 Sco Kbl D l (1 + z)~ 3 > 

where L' co is the luminosity, v b s is the observed 
frequency in units of GHz, Dl is the luminosity 
distance in Mpc and z is redshift. Finally molec- 
ular hydrogen mass Mr 3 is calculated as: 

Mh 2 = L' co x a C o, 



where aco is the Galactic conversion factor of 
3.2M Q \(K km sec~ 1 pc 2 )~ 1 ] which does not in- 
clude a correction for the presence of Helium. 

3. Results 

Combining the data from the diverse set of ob- 
servations presented in Section 2, the following 
picture on the nature of VGS_31 emerges. 

VGS_31 consists of three galaxies located in a 
void, VGS_31a, b and c (Figure 1). A careful study 
of the HI properties suggests that the galaxies are 
embedded in an elongated HI cloud (Figure 3 and 
4). In addition, both the optical and the HI obser- 
vations suggest strong interactions in the system 
(Figure 3). This emerging picture tells us that we 
may be dealing with two different processes. One 
is the assembly of a filamentary structure in a void, 
the other is an interaction between the galaxies. 
We will describe these two main results separately. 

VGS_31 exhibits a very peculiar structure in 
HI. The HI column density map and the position- 
velocity (PV) diagram (Figure 3 and 4) present 
two views elucidating the nature of the system; 
i) the filamentary structure and ii) its kinematic 
properties. 

3.1. Filamentary structure 

The PV diagram seen on the right hand side 
of Figure 4, is created by taking one single slice 
through the HI data cube, from the far east to the 
far west of the VGS_31. In the HI column density 
map (Figure 3) galaxies appear to be lined up in a 
HI filament. The whole HI envelope extends over 
~ 120 kpc on the sky. The PV diagram shows that 
the three galaxies are at about the same systemic 
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Fig. 7.— IRAM spectrum of VGS_31b in the 
115.271 GHz. The intensity scale is in T%_, in 
Kelvin. 
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Fig. 8. — H Q major and minor axis spectra and images for VGS_31b (left) and VGS_31a (right). In the two 
figures above, the H Q emission line spectra along the optical major and minor axis of VGS_31b and VGS_31a 
are shown on the same scale as their optical images. In both figures, the optical images of the galaxies are 
overlaid by white H a imaging contours. On top of these, the yellow lines indicate the positions of the slits. 
b represents the slit position along the optical major axis and a represents the slit position along the optical 
minor axis. The optical images of the galaxies and the spectra are scaled in such a way that the physical 
sizes of the line profiles matches with the real size of the corresponding H Q emission region within ~ 0.6 kpc. 
Each diagram showing optical spectra have axes of relative distance from the centre (in kpc) vs. heliocentric 
radial velocity (in km/s). 



Table 3 
h q fluxes and luminosities 



Galaxy 


hg[F(H a + Nil)] 


\og[[F(H a )} 


log[I(H a )} 


L(H a ) 
{erg s _1 ) 


VGS_31a 


-12.622 


-12.755 


-12.304t0.076 


4.7 xlO 41 


VGS_31b 


-12.345 


-12.478 


-11.913t0.104 


1.16xl0 42 


VGS_31c 


-13.604 


-13.667 


-13.527t0.019 


2.8 xlO 40 



Note. — Column 1: Galaxy name. Column 2: Measured H Q flux including 
Nil lines. Column 3: H a flux corrected for Nil deblending. Column 4: Extinc- 
tion corrected H a flux. Column 5: H a luminosity derived from the extinction 
corrected H Q flux. 
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Table 4 
Star formation properties. 



Galaxy 


SFR a 
{M Q yr- 1 ) 


SFRjjy 

(Moyr- 1 ) 


S S F Rq, 
(lO-^oyr- 1 ) 


SFR a /M HI 
(lO- s M yr- r ) 


SFR a /M H2 
(10-SMoyr- 1 ) 


VGS_31a 


2.5 


1.4 


0.07 


0.13 




VGS_31b 


6.3 


2.2 


0.06 


0.43 


0.7 


VGS_31c 


0.15 


0.22 


0.05 


0.09 





Note. — Column 1: Galaxy names. Column 2: Star formation measured from H a flux lumi- 
nosities. Column 3: Star formation calculated from NUV Column 4: Specific star formation 
rate. Column 5: Star formation rate per HI mass. All three galaxies have similar SSFRaS 
while their SFR a /Mni are different. 



Table 5 

Stellar, HI and molecular masses and CO(1-0) luminosities for VGS31. 



Galaxy 


M, 


M H i 


Leo 


M H2 


M H JM, M h JM hi 




(1O 8 M ) 


(1O 8 M ) 


(10 8 K km s-V 2 ) 


(1O 8 M ) 




VGS_31a 


35.1 


19.89 ± 2.90 








VGS_31b 


105.31 


14.63 ± 1.97 


2.8 


8.9 


0.085 0.608 


VGS_31c 


2.92 


1.66 ± 0.95 









Note. — Stellar masses are taken from the publicly available MPA-JHU catalog for SDSS DR7 
which were derived using fits to the broad-band ugriz photometry. VGS31 galaxies are all spectroscopic 
targets in SDSS DR7 with 3" fiber spectra in their central regions and these spectra have been analyzed 
and included in this catalog. 
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velocity. Also the extreme ends of the hiament are 
at the same velocities suggesting a single, coherent 
filamentary structure. 

3.2. Kinematic properties 

The kinematics along this perceived filamentary 
structure is, however, rather complex. We will 
examine this, galaxy by galaxy, going from east to 
west in Figure 4. 

VGS-31b: The optical tail has a HI counter- 
part, showing almost no velocity gradient. There 
is a slight offset between this HI tail and its optical 
counterpart. It is kinematically connected to the 
inner disk and the ring of VGS_31b. The spatial 
resolution is insufficient to distinguish the gas as- 
sociated with the ring from the gas associated with 
the disk. Therefore, we can not say whether the 
HI tail is connected to the inner disk or the ring. 
The galaxy disk and the ring exhibit rotation with 
a velocity spread of 325 km s _1 . The H a spectrum 
(Figure 8) shows the rotation along the optical ma- 
jor axis. In the PV diagram, the velocity gradi- 
ent is not as steep because the slice has not been 
exactly placed along the optical major axis but 
rather towards the tail. West of the disk, we see a 
HI bridge between VGS_31b and VGS_31a. Close 
to VGS_31b the HI in this bridge has a velocity 
dispersion of ~ 50 km s _1 . In VGSSlb, there 
is almost no velocity gradient in the H a emission 
along the optical minor axis whereas we see a sig- 
nificant gradient along the optical major axis. 

VGSSla: However, halfway to VGS_31a the 
velocity dispersion becomes significantly larger, ~ 
125 km s . This behaviour is repeated at the 
other side of VGS_31a as well. The HI bridge 
has not been detected in the optical or UV. The 
gas associated with VGS_31a shows broad veloc- 
ity range. However, we don't have the spatial 
resolution to map the HI kinematics in detail. 
As we move away from the galaxy towards the 
west, we continue to see HI with a broad veloc- 
ity dispersion. The H Q spectrum (Figure 8), has 
much higher resolution and shows a large kine- 
matic asymmetry. 

VGSJSlc: Very little HI is detected between 
VGS_31a and VGS_31c. A clear HI connection is 
however visible in Figure 4. The velocity disper- 
sion of the HI bridge becomes again narrower as 
we approach VGS_31c. As for the bridge between 



VGS_31b and VGS_31a, there is no stellar counter- 
part for this HI bridge either. The HI resolution 
is not enough to determine the detailed gas kine- 
matics of VGS_31c. 

3.3. Structure of the galaxies 

In addition to HI, optical data complement the 
picture of VGS_31's complex dynamics. We will 
continue to describe the peculiar features of the 
system going from east to west. 

VGSJSlb: As seen from Figure 1, VGS_31b is, 
optically, the most disturbed galaxy. It has a one 
sided tail, curved towards the north east. There is 
a ring like structure around the disk. The tail and 
the disk seem to be connected. The first zoomed 
image shows the inner disk. The second zoomed 
image displays the bar positioned asymmetrically 
in the disk and the bright central part. The H Q 
emission of VGS_31b shows clear rotation along 
the optical major axis. In addition there is a very 
fast rotating inner structure, about 1 kpc in ex- 
tent with a velocity width of ~ 500 km/s. Along 
the minor axis a minor velocity diference is seen, 
indicative of non-circular motions along the bar 
(Figure 8). Also, all star formation is concentrated 
in this central region where the bar is present. 
VGS_31b has been detected in CO(l-O). The ve- 
locity of the CO (1-0) peak is at ~ 6200 km s _1 , 
which is similar to the velocity of the HI. 

VGSSla: The zoom-in image (Figures 1) of 
VGS_31a shows that it is slightly disturbed and 
the bar like structure in the center overlaps with 
the H Q emission region (Figures 5). Like VGS_31b, 
rotation is along the optical major axis. However, 
the H Q emission line profile along the major axis 
is quite irregular. 

VGS-31c: This is the smallest of the three 
galaxies without any significant morphological ir- 
regularities (Table 5). As in the other two galax- 
ies, the star formation is confined to the central 
part of the disk. 

3.4. Star formation properties 

The H a and UV results show that all three 
galaxies exhibit recent star formation activity con- 
centrated in their central parts (Figure 5 and Ta- 
ble 4 ). The tail and the ring of VGS_31b are 
not detected in H a or UV. The same is true for 
the HI bridges between VGS_31b, VGS_31a and 
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VGS_31c. Among the three galaxies, VGS_31b 
has significantly higher SFR Q (Table 4). We also 
emphasize that the small D„(4000) break values 
are indications for young stellar populations in the 
central part of the galaxies. 

From the SDSS spectra we may determine the 
location of the VGS_31 members in a Baldwin, 
Phillips & Terlevich (BPT) diagram (Baldwin 
et al. 1981) shown in Figure 6. For this diagram, 
emission lines of all galaxies, including VGS_31, 
have been extracted from the SDSS DR7 spectral 
database of 3" fiber apertures. In this diagram, 
both VGS_31a and VGS_31b are located in the 
HII zone together with the other starburst galaxy 
samples. Their SDSS spectra have strong Balmer 
emission lines with large equivalent widths and 
blue continua. VGS_31a and VGS_31b are located 
inside the HII/ starburst region in the BPT dia- 
gram while VGS_31c is placed inside the normal 
star forming zone. 

Stellar masses of the three galaxies range from 
3 x 1O 8 M to 1.06 x 1O 1O M (Table 5). In con- 
trast to the difference in their stellar masses, they 
have similar S_SFR a . It is worth nothing that 
their S_SFR Q and SFR q /Mhi are significantly 
above the median of those of the ALFALFA aver- 
age density sample (Kreckel et al. 2012), indicat- 
ing enhanced star formation in comparison with 
galaxies of similar mass and gas content. 

4. Discussion 

VGS_31 is a peculiar system through which we 
may be witnessing the ongoing growth of three 
galaxies along a filament inside a void. Here we 
discuss possible scenarios for the evolution of the 
VGS_31 galaxies in the void including gas accre- 
tion from an intra- void filament, interactions and 
merging. These processes are difficult to disentan- 
gle from one another on the basis of the observa- 
tional material presented. Yet we will discuss the 
observed phenomena in the context of all of these 
processes and indicate which we consider most im- 
portant in each of the individual galaxies. 

First we will discuss different interaction sce- 
narios for VGS_31a and VGS_31b using the ob- 
servational results presented above. VGS_31b is 
the most eye-catching member in the system. A 
one sided tidal tail and a ring like structure are 
clearly visible in the optical (Figure 1). This mor- 



phology suggests a minor merging incident with 
a low mass galaxy. A satellite galaxy could have 
left the tail and formed the ring by wrapping the 
disk of VGS_31b (Mihos & Hernquist 1994, 1996; 
Due & Renaud 2011). It is less likely that the 
tail has been caused by tidal interaction between 
VGS_31b and VGS_31a. As the mass ratio of these 
systems is 3 to 1 one would expect greater damage 
to the disks and more prominent tails and counter- 
tail features as usually seen in tidal interactions 
and major mergers (Hibbard & van Gorkom 1996). 
The disks of VGS_31b and VGS_31a are not de- 
stroyed as expected in minor mergers (Schweizer 
2000; Martmez-Delgado et al. 2010; Due & Re- 
naud 2011). 

VGS_31b has a bar, visible in close up images 
in Figure 1. The bar is pronounced in H a as well 
(Figure 5). The kinematics of the H Q (Figure 8) 
indicate a fast rotating inner structure and evi- 
dence for streaming motions characteristic for a 
bar. In addition, VGS_31b is a starburst Markar- 
ian galaxy having enhanced star formation in its 
central part overlapping with the bar (Figure 5 
and 8, Table 4). The detection of CO and loca- 
tion of VGS_31b in the BPT diagram support the 
starburst picture, presumably the consequence of 
gas accretion into its central part as seen in many 
galaxy mergers and interactions (Mihos & Hern- 
quist 1994, 1996; Due & Renaud 2011). There is 
an offset between the stellar and the gas compo- 
nent of the tail (Figure 3). This is also observed 
in some tidally interacting systems (Mihos 2000) 
and may be due to different initial distributions 
of both components or additional processes that 
act on one component and not on the other (see 
Due & Renaud 2011, for a recent review). On the 
other hand, in most of the tidal tails in mergers, 
ongoing star formation is observed (Hibbard & 
van Gorkom 1996; Ncff et al. 2004), in the case 
of VGS_31b, however, no ongoing star formation 
has been detected neither in the tail nor in the 
ring (Figure 5). 

VGS_31a, however, shows different characteris- 
tics. It has no visible tails, but the internal kine- 
matics is disturbed (see the H Q kinematics in Fig- 
ure 7) and the HI shows that gas east and west 
of the galaxy exhibits a large spread in velocity as 
if there is a corotating halo filled with HI rather 
than tidal features with simple kinematics (Fig- 
ure 4). Like VGS_31b, VGS_31a is a starburst 
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galaxy with enhanced star formation in its central 
part (Figure 5 and 6, Table 4). There are two pos- 
sibilities for the mechanism which could cause the 
enhanced star formation and the morphological 
disturbance: i) VGS_31a could be in interaction 
with VGS_31b and this may cause the disturbence 
in its morphology and results in the irregular H a 
emission line profile, ii) Instead of accreting gas 
from VGS_31b, VGS_31a may experience steady 
gas infall from the intergalactic medium, presum- 
ably from the structure outlined by the location 
of the galaxies and the HI 'connecting' them (Fig- 
ure 4). This could cause the enhanced star for- 
mation as in the first scenario and better explains 
the broad velocity range of the HI surrounding 
VGS_31a. One could argue that if VGS_31a ex- 
hibits such accretion from an LSS filament, then 
one would also expect VGS_31b to show the same 
characteristics. On the other hand it is conceiv- 
able that a minor merger, as suggested by the tail 
and ring could have disrupted the accretion pro- 
cess. Detailed simulations of such scenarios are 
required to test the validity of these scenarios. 

VGS_31c, the smalles of the three galaxies, has 
enhanced star formation and its HI shows the 
charecteristics of interactions, albeit at very low 
signal to noise. It is difficult in this case to deter- 
mine conclusively the process (es) responsible for 
the HI structure and enhanced star formation. 

The most exciting result is that we may be wit- 
nessing the assembly of structure within a void, 
and the birth process of galaxies in such a deso- 
late area. It does fit into the theoretically expected 
buildup of voids and galaxies therein. Voids evolve 
in a hierarchical fashion, leaving planar and fila- 
mentary substructure within the emerging voids 
(Dubinski et al. 1993; van de Weygaert & van 
Kampen 1993; Sahni et al. 1994; Sheth & van 
de Weygaert 2004; Aragon-Calvo & Szalay 2013) 
and (Rieder et al. 2013 (subm.)). The question 
is whether we see here the manifestation of this 
process. VGS_31 could be a density enhancement 
within an underlying dark matter filament. The 
complication is that in addition to accretion of ma- 
terial from a LSS filament the galaxies also suffer 
from tidal interactions and minor merging. As for 
the cold flow accretion scenarios we can only make 
a very rough estimate because we don't know ex- 
actly the orientation of the system. If we take the 
size of the HI filament and the velocity gradient 



from east to west then we can estimate a rough 
accretion timescale of at least a gigayear. An im- 
portant next step will be to use advanced simula- 
tions with gas and star formation to see whether 
the scenario proposed here does indeed take place 
in voids. 

Finally, it is of interest to note that VGS_31 
is not the only case for a filamentary structure 
and cold flow accretion in our Void Galaxy Survey. 
We have discovered at least two other cases which 
seem to exhibit this hierarchical structure forma- 
tion. In addition to another filamentary galaxy 
configuration, VGS_38 (Kreckcl et al. 2011), we 
also found a polar disk galaxy, VGS_12 (Stanonik 
et al. 2009). VGS_38 is a system of chain galaxies 
which share the same HI. VGST2 is located right 
at the center of a tenuous wall between two large 
voids. It has a polar HI disk much more extended 
than its stellar disk. The polar disk has no stellar 
counterpart or any ongoing star formation. This 
galaxy is a candidate for the cold flow accretion. 

Individually these galaxies and VGS_31 are un- 
usual systems, however taken collectively they 
show that the void environment is an extremely 
interesting site for understanding galaxy forma- 
tion and evolution. 
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